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The infinite disorder fixed point of the random transverse-field Ising model is expected to control 
the critical behavior of a large class of random quantum and stochastic systems having an order 
parameter with discrete symmetry. Here we study the model on the square lattice with a very 
efficient numerical implementation of the strong disorder renormalization group method, which 
makes us possible to treat finite samples of linear size up to L = 2048. We have calculated sample 
dependent pseudo-critical points and studied their distribution, which is found to be characterized 
by the same shift and width exponent: v = 1.24(2). For different types of disorder the infinite 
disorder fixed point is shown to be characterized by the same set of critical exponents, for which we 
have obtained improved estimates: x = 0.982(15) and ?/) = 0.48(2). We have also studied the scaling 
behavior of the magnetization in the vicinity of the critical point as well as dynamical scaling in the 
ordered and disordered Griffiths phases. 



I. INTRODUCTION 



The random transversc-ficld Ising (RTFI) model is the 
prototype of random quantum systemai having a quan- 
tum critical point at zero temperature^. This model has 
experimental realizations^ and there is a large amount of 
theoretical work, which aims to clarify the properties of 
the random critical point. It is expected that basic fea- 
tures of the critical behavior are demonstrated in the onc- 
dimcnsional [Id) model and therefore most of the theo- 
retical studies are performed in Id. After early works by 
McCoy and others^ Fisher— has used a renormalization 
group (RG) framework to obtain several presumably ex- 
act results about its critical properties. It has been shown 
that the critical properties of the Id model are governed 
by a so called infinite disorder fixed point (IDFP), in 
which the strength of disorder grows without limit dur- 
ing renormalization^. As a consequence disorder fluctu- 
ations are dominated over quantum fluctuations and the 
approximations used in the RG approach become exact at 
the critical point. The IDFP of the RTFI model is shown 
to govern the critical properties of another random quan- 
tum systems having an order parameter with discrete 
symmetry^i^. Furthermore this fixed point is found to 
be isomorphic with that of several Id stochastic models, 
such as the Sinai walki^, the random contact proccssii 
or the random exclusion process^^. 

Comparatively less results are available about the crit- 
ical behavior of the RTFI model in higher dimensions, 
which arc almost exclusively restricted to 2d. By now 
different numerical studies are in favor of the conclu- 
sion, that also in 2d the critical behavior is controlled 
by an IDFP. In this respect we mention different numeri- 
cal implementations^^"— of the strong disorder renormal- 
ization group (SDRG) method, as well quantum Monte 
Carlo simulationsi^. These results are in agreement 



with recent simulation studies of the 2d random contact 
process^^, which is expected to belong to the same univer- 
sality class. Also the 2d random walk in a self-affine ran- 
dom potential^! could be related to the 2d RTFI model. 
As far as the numerical estimates of the critical exponents 
in 2d are concerned, these contain quite large errors, for 
a summary of the estimates see Ref.fl^l- In simulation 
studies these are connected to the logarithmically slow 
critical relaxation, whereas in the SDRG method the er- 
rors has mainly finite-size origin the typical linear size of 
the largest systems being about L = 128— 160. Also the 
type of disorder used in the calculations has an influence 
on the error of the results. Due to these numerical limi- 
tations there are no quantitative results in 3d, although 
it is very probable that the random critical point is an 
IDFP in this case, too^^. 

In this paper we are going to revisit the problem of 
the critical behavior of the 2d RTFI model. Like other 
studies we use numerical implementation of the SDRG 
method, however we have developed a very efficient algo- 
rithm, which make us possible to treat systems as large 
as L = 2048. In this way the number of sites in our sam- 
ples are several hundred larger, than in previous studies. 
Comparing with earlier SDRG investigations our study 
has several different features, i) We define and calculate 
finite-size pseudo-critical points and study their distribu- 
tion, ii) We obtain accurate estimates for the true criti- 
cal point of the model, calculate effective, size-dependent 
critical exponents and study their extrapolation, iii) We 
consider different forms of the disorder and study the uni- 
versality of the critical exponents as well as the scaling 
functions, iv) We also study scaling outside the critical 
point, as well as dynamical scaling in the disordered and 
ordered Griffiths phases. 

The structure of the rest of the paper is the following. 
The model and the SDRG method is presented in Sec. 
im The basic features of the computer algorithm arc dis- 



cussed in Sec. IIIIl In Sec. IIVI we describe how finite-size 
transition points are defined and calculated within the 
frame of the SDRG method. We study their distribution 
and analyze the shift of the mean value as well as the 
width as a function of the size of the system. In Sec. |V] 
critical exponents are extracted through finite-size scal- 
ing and the scaling behavior outside the critical point is 
analyzed. We also study dynamical scaling in the disor- 
dered and in the ordered Griffiths phases. Our results 
are discussed in the final section, whereas some details of 
the computer algorithm are presented in the Appendix. 



II. MODEL AND THE SDRG METHOD 

A. Random transverse-field Ising model 

We consider the RTFI model defined by the Hamilto- 
nian: 
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J2 ^^^'^ 



(1) 



in terms of the Pauli- matrices, fr^^'^. Here i and j are sites 
of a square lattice and the first sum runs over nearest 
neighbors. The Jij couplings and the hi transverse fields 
are independent random numbers, which are taken from 
the distributions, p{J) and q{h), respectively. 

Here we use two different type of distributions, which 
have both the same uniform distribution of the couplings: 



p(j) = e(j)e(i-j). 



(2) 



Q{x) being the Heaviside step-function. For the 'box-/i' 
disorder we have: 



qih) 
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e{h)e{hb-h) 



(3) 



whereas for the 'fixed-ft,' model we have a constant trans- 
verse field: 



q{h) = S{hf - h) 



(4) 



In the following we use the logarithmic transverse-field 
= In /if, or = In hf to characterize the system. In 
the thermodynamic limit, L —> oo, the system in Eq.((T|) 
displays a paramagnetic phase, for 9 > 9c, and a ferro- 
magnetic phase, for 9 < 9c- In between there is a random 
quantum critical point at 9 ~ 9c- The quantum control- 
parameter is defined as 5 = 9 — 9c- 



and at the same time new terms are generated between 
remaining sites by second-order perturbation method. 
The procedure is sketched in Fig. [T] for a higher dimen- 
sional system. If the largest term is a coupling (see the 
right panel of Fig. [T]), say il = Jij connecting sites i and 
j, then the two sites involved form a spin cluster with an 
effective moment fi' — fii + fij (initially /i^ = 1 Vi) in the 
presence of an effective transverse field: h' w hihj/Jij. 
The renormalized value of the coupling of the cluster to a 
site a is given by the 'maximum rule' max [Ja,i, Ja,j]- On 
the other hand, if the largest local term is a transverse- 
field (see the left panel of Fig. [1]) , say fl = hi, then site 
i is eliminated and effective couplings are generated be- 
tween each pairs of spins, among the neighbors of i. If 
a and b arc neighboring spins to i, then the generated 
coupling is given by: J'^^ i, ~ Ja,iJb.i/hi- If the sites a and 
b are already connected by a coupling, Ja.t 7^ 0, than the 
renormalized coupling is given by the 'maximum rule' as 
max [Ja,b, J'a b\- T^^G use of the maximum rule is justified 
if the renormalized couplings have a very broad distri- 
bution, which is indeed the case at the IDFP. We shall 
see that with the maximum rule the numerical algorithm 
can be simplified. At each step of the renormalization 
one site is eliminated and the energy scale is continu- 
ously lowered. For a finite system the renormalization 
is stopped at the last site, where we keep the energy- 
scale and the total moment as well as the structure of 
the clusters. 



ft= J; 





J' — JaiJbi/hi 

J = ma.x{Jad, JaiJdi/hi) 



J — max(jQj, Jaj) 



FIG. 1: Illustration of the decimation steps of the strong 
disorder renormalization group method in higher dimensions. 



B. Strong disorder renormalization group metiiod 

Here we use the SDRG method^^, which has been in- 
troduced by Ma, Dasgupta and Hu^ and later developed 
by D. Fisher— and others. During the method, the largest 
local term in the Hamiltonian (which defines the energy- 
scale, ri, at the given RG step) is successively eliminated 



C. Scaling at the infinite disorder fixed point 

At the IDFP the distribution of the effective couplings 
and that of the transverse fields becomes broader and 
broader during the renormalizatior*2i^. Considering the 
ratio of two effective terms in the Hamiltonian at a given 
stage of the SDRG it will tend either to infinity or to 
zero. This indicates that the disorder is infinitely strong 



and the pcrturbative results during the RG are exact. 
Quahtatively the log-energy scale, lnJ7, scales with the 
linear size of the system size, L: 



\n{no/n) r^L"^ 



(5) 



where Oo denotes a reference energy scale. The aver- 
age spin-spin correlation function is defined as G{r) = 
{<7i<^i+r)^ where (...) denotes the ground-state average 
and (...) stands for the averaging over quenched disor- 
der. The asymptotic value of the correlation function 
defines the magnetization, m, in the system: 



lim G{r) = m 



(6) 



and TO > in the ferromagnetic phase and tti = in the 
paramagnetic phase. The connected correlation function, 
G{r) = G{r) — m? , in the vicinity of the critical point 
behaves as: 



G{r) - r^2:r exp(-r/^) 



(7) 



where the correlation length, ^, is divergent at the critical 
point as: 



i-\5[ 



(8) 



Thus at (5 = there is a power-law decay of the corre- 
lations, which is related to the fractal structure of the 
spin clusters. Indeed, the average cluster moment, /i, is 
related to the energy-scale, £1 as: 

^l ^ [Hno/n)X^ , (9) 

and can be expressed also with the size: 

^i^L'^f . (10) 

Here the fractal dimension of the cluster, df, is related 
to the other exponents as: 

df = 0V = d — X . (11) 

In Id the critical exponents are exactly knowiiii^: 



1 



-^^, ^ = 2, 2- = ^-. (12) 



during which numerous new bonds arc generated and as 
a result sites with large number of links are formed. In 
this way, after a nai've implementation of the method, 
the system is transformed into an almost complete graph 
and the subsequent ^.-decimations, generating new links 
between practically all pair of sites, will be very slow. For 
a system with N sites this algorithm would need 0{N^) 
time. 

Using the maximum rule in the approach, however, of- 
fers two ways to speed up the procedure^. First, one can 
notice that the renormalization trajectory is not unique 
in this case. There are terms in the Hamiltonian, which 
are called as "local maxima" and which can be decimated 
independently. Thus one should not follow the "decima- 
tion of the largest term in each step" principle, instead we 
are going to optimize the time of the renormalization tra- 
jectory, which goes over in some order of the local max- 
ima. The second consequence of the maximum rule is, 
that a large number of bonds will never be participating 
in the renormalization process. These latent bonds can 
be deleted from the list of edges without consequences. 
The latent bonds are in such a local environment, that 
after decimating a nearby site or bond a stronger new 
coupling is generated to the same edge, thus the original 
bond disappears without participating in the renormal- 
ization. Filtering out these irrelevant bonds will result in 
a considerable improvement of the algorithm. In the fol- 
lowing we discuss the properties of the local maxima and 
the optimal RG trajectory as well as the main features 
of the filtering process. 



B. Local maxima and the optimal RG trajectory 

A local maximum in the set of couplings and transverse 
fields is such a term, which is larger (not smaller) than 
any of its neighboring terms. Considering a coupling Jy 
is a local maximum, provided Jij > hi, Jij > hj, and 
Jij > Jiki V/c, as well as Jij > Jij, VL Similarly a trans- 
verse field, hi, is a local maximum, if hi > Jij, Vj. It 
can be shown that the local maxima can be decimated 
independently, the renormalization performed in any se- 
quence gives the same final result. 



C. Filtering out irrelevant bonds 



III. IMPLEMENTATION OF THE SDRG 
METHOD 

A. Problems with the naive implementation 

The SDRG decimation rules are very simple and it is 
straightforward to implement the method numerically. In 
higher dimensions, however, the topology of the lattice 
is changing during renormalization, which could result in 
considerable increase of the computational time. More 
dangerous steps in this respect are the /i-decimations. 



The principle, that some bonds are irrelevant and 
does not modify the renormalization procedure has been 
first realized by Kawashima— . He has also suggested 
a criterion to identify the irrelevant bonds which are 
then deleted from the graph. This filtering procedure 
as illustrated in Fig. [5] is used in a few 2d numerical 
^Qj-]jgi5ji6ji8 jjj Appendix [^ we give the proper defini- 
tion of the filtering criterion and prove it. 

The use of the filtering criterion for a site with k neigh- 
bors needs typically 0{k^) time, since 0{k^) triangles 
have to be checked. However, the application of the fil- 



Before decimation 



After decimation 



After filtering 






FIG. 2: Illustration of the filtering criterion. 



suming. In step 3. making the filtering before decima- 
tion (prefiltering) will ensure that the dangerous bonds 
are deleted. 

This algorithm in 2d works typically in 0{L^ In L) time 
in a i X L system near the critical point. With our 
method an L = 128 sample is renormalized in ^ 1 second, 
whereas for L = 1024 the typical time is ^ 1.5 minutes 
(in a 2.4GHz processor). 



tering makes not too much reduction of the computa- 
tional time, if the bond is in a general position, where 
not too many new terms arc generated after consecutive 
RG steps. There are, however, bonds in "dangerous posi- 
tions" , having at least one of the neighboring transverse 
fields as a local maximum. Decimating this local max- 
imum many new terms are generated and therefore the 
use of prefiltering, which checks the bonds of the deci- 
mated spin just before its decimation is very effective, 
which typically needs only 0{k^) time. 



D. Basic steps in the numerical SDRG method 

In order to obtain an efficient implementation of the 
SDRG method we combine the optimal selection of the 
RG trajectory with the filtering algorithm. In our 
method we use the following steps. 

1. We check all terms of the Hamiltonian and select 
the local maxima. We make two lists, one for the 
couplings and one for the transverse fields. 

2. We decimate every coupling, which is on the list of 
local maxima. In this step the order of the decima- 
tion is arbitrary. After decimating the original cou- 
plings new terms arc generated, among which there 
are new local maxima. We include those into the 
list of local maxima and at the same time we dec- 
imate the new couplings from this list. We repeat 
this step, until the list of local maximum couplings 
is empty. 

3. If all or all but one the transverse fields are local 
maxima we select the smallest one, delete all the 
rest sites and END the iteration. Otherwise we se- 
lect a transverse field from the list of local maxima 
having a small (or the smallest) degree, which is de- 
fined as the number of edges to the given site and 
filter for the neighboring bonds. We decimate this 
transverse field and check the generated new terms 
for local maxima. If the list of local maximum cou- 
plings is not empty, we go to step 2. Otherwise we 
go to step 3. 

Using the selection rule in the second part of step 3. will 
ensure that the average degree of the sites will not be 
large. In this way we prevent the formation of too con- 
nected clusters, the decimation of which being time con- 



IV. FINITE-SIZE CRITICAL POINTS 

A. Scaling of pseudo-critical points 

In the study of the critical behavior of random systems 
it is very important to find an accurate estimate of the lo- 
cation of the critical point. The quality of the estimate of 
9c will influence the error of the calculated critical expo- 
nents and scaling functions. In a random sample of linear 
size, L, one can generally define finite-size pseudo-critical 
points^^"— , Oc{L), which are usually given as the position 
of the maximum of some physical quantity, which is di- 
vergent in the thermodynamic limit (c.f. susceptibility) 
at 6c- The distribution of 6c{L) provides important in- 
formation about the scaling behavior at the fixed point 
of the system3^. In particular one studies the shift of the 
average value, Oc{L), which is expected to scale as: 



h{L) 



L 



^Iju, 



(13) 



with the shift exponent, Vs- Similarly, one measures 
the width of the distribution, AddL), which behaves for 
large-L as: 



AdciL) ^ L- 



^i/i^^ 



(14) 



where z/^, denotes the width exponent. According to 
renormalization group thcory^^ for a classical random 
system with relevant disorder— the critical behavior is 
controlled by a conventional random fixed point, with 
i^s ~ Vw ~ V, where v > 2/S^ is the correlation length 
critical exponent of the system. 



B. Identification of pseudo-critical points 

In a random quantum system one has to use another 
methods to locate pseudo-critical points^. One method, 
which is well suited to the SDRG approach is the dou- 
bling method^, which has been used for chains^ and 
for laddersi^ of the RTFI model. In the doubling pro- 
cedure in Id (or quasi- id) geometry one considers a 
random sample (a) of length, L, and makes a dupli- 
cated sample (2a) of length 2L by joining two copies 
of (a). Using the SDRG method one calculates some 
physical quantity (magnetization or gap) in the origi- 
nal and in the replicated sample, which is denoted by 
f{a,L) and f{2a,2L), respectively, and study their ra- 
tio, r{a,L) = /(2a, 2i)//(a, L), as a function of the 



control parameter, 9. At — 9c{a, L) this ratio has a 
sudden jump, which is identified with the pseudo-critical 
point of the sample. We note that the actual value of 
Oc{a, L) is practically independent of the physical quan- 
tity we consideredi^, since this singularity is connected 
to the topology of clusters produced by the SDRG ap- 
proach. 

Having this observation in mind we can generalize the 
doubling method for two (and higher) dimensions. In 2d 
we glue together two identical square-shaped samples at 
the boundaries as indicated in Fig. [3] Then we renor- 
malize the duplicated sample and calculate the structure 
of the connected clusters, among which there might be 
such, which have sites (in equivalent positions) in both 
replicas. These sites are correlated and the fraction of 
these correlated sites defines the correlation function be- 
tween the replicas. By increasing the control parameter, 
9, the replica correlation function is decreasing and at a 
well defined value, 9c{L), it suddenly jumps to zero. We 
consider 9c{L) as the pseudo-critical point of the given 
sample. It is easy to see that this definition is equiva- 
lent to the previously used criterion in Id, furthermore 
it is straightforward to generalize it to three- and higher 
dimensions. 




FIG. 3: Illustration of the boundary conditions used in the 
doubhng procedure in higher dimensions. 



C. Numerical results 

1. Distribution of pseudo-critical points 

We have calculated pseudo-critical points of square- 
shaped samples for various linear sizes, L, which are ex- 
pressed as L = 2' and L = 3 * 2'^^ up to / = 10. Gener- 
ally we have considered 4 x 10"* realizations except for the 
largest system, when we had at least 10^ samples. The 
distribution of the pseudo-critical points is shown in Fig|3] 
for both type of disorder. The mean value of the critical 
points is considerably larger for box-ft, randomness and 
also the width of the distribution - for the same value of 
L - is larger in this case. (We note that the same trend 
is seen in Id, where 0c = and 9c = —1.) Taking 
into account the result in Ea.(|13p the appropriate scal- 
ing combination is y = {9c{L) — 9c)L^''^ in terms of which 
the scaled distributions, p{y), are shown in the insets of 
FigHl Here using our final estimates in Egs. ip^ and (PT|) 



we obtain excellent scaling collapse of the data points for 
both type of randomness. The scaling curves for the two 
different disorder approach the same standardized mas- 
ter curve, which indicates that the fixed point of the RG 
transformation is unique and (at least for strong enough 
disorder) strongly attractive. The master curve is differ- 
ent from that in Id, which is Gaussian in this caso^. In 
2d the maximum of the curve is shifted to negative values 
and the distribution is non-symmetric. We have calcu- 
lated the percolation (or spanning) probability, Ppr, at 
the critical point, which is given by the fraction of sam- 
ples having finite replica correlation function at 9c. It 
can be expressed with the scaled distribution function as 
Ppr — /g p{y)dy. Our estimate is: 



Ppr = 0.149(2) 



(15) 



which is much smaller, than for standard 2d 
percolation^. We have also calculated the skewness, s, of 
the distribution, which has the value s = 0.19(3) for both 
type of disorder. The asymmetric form of the distribu- 
tion indicates that in 2d the topology of the renormalized 
model is different from that in Id. Samples having more 
strongly connected clusters, thus a higher 9c{L), have a 
somewhat larger weight than the less strongly connected 
clusters. 



2. Shift of the finite-size critical points 

For a fixed linear size, L, we have calculated the mean 
value, 9c{L), as well as the median 9^°'^{L) of the pseudo- 
critical points. Since the distribution is non-symmetric 
these two characteristic values are not identical, however 
both are expected to follow the scaling form in Eq.(fT3|) 
with the same value of the shift exponent, t^^. This is 
illustrated in FiglSl in which 9c — 9c{L) as well as 9c — 
^mcd^^j is shown as a function of L in a log-log scale for 
both type of disorder. Having appropriate limiting values 
for 9c, see Eg. ip^ . the points in Fig jS] are asymptotically 
very well on straight lines (both for the mean value and 
for the median) having approximately the same slopes: 
-1/iy, « -0.8. 

To get more quantitative estimates we have calculated 
effective, size dependent shift exponents which are de- 
fined as: 



i^siL) 



ln2 



■In 



:(2L) - 9ciL) 



,{L)-9c{L/2) 



(16) 



and similarly for the median values. These are shown in 
the two insets of Fig 15] as a function of 1/L (upper in- 
set for the mean) and InL (lower inset for the median), 
respectively. The exponents calculated from the mean 
values show 1/L correction terms for both type of dis- 
order, although with different signs. The extrapolated 
values are 1/i/f' = 0.79(2) and l/i^P == 0.81(2), which 
agree within the error of the calculation and leads the 
estimate: l/iyg — 0.80(2). As seen in the lower inset of 
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FIG. 4: (Color online) Distribution of the pseudo-critical 
points, 0c{L), for various lengths for hox-h randomness (up- 
per panel) and for fixed- /i randomness (lower panel). In the 
insets the scaled distributions are shown as a function of 
y = {OciL) - ec)L^^", see the text. 



Fig IS] the effective exponents calculated from the median 
of the distribution have weaker 1/L dependence, instead 
they show log-pcriodic-likc variations. The estimates for 
l/vs from these data are compatible with the previous 
estimates obtained from the mean values, thus wc can 
write our estimate of the shift exponent: 



,^s = 1.25(3) 



(17) 



3. Scaling of the 



wi 



idth of the pseudo-critical points 



We have measured the standard deviation of the dis- 
tribution of the pseudo-critical points, A9c{L), which are 
shown in Fig. |6] as a function of the linear size, L, in 
a log-log scale. In agreement with the scaling predic- 
tion in Ea. ((Ti|) the points in FiglHlare asymptotically on 
straight lines, the slope of which is approximately the 




FIG. 5: (Color online) Scaling of the shift of the finite-size 
transition points, calculated from the mean value (+ for 'b' 
and ^ for 'f'), as well as from the median ( for 'b' and 
El for 'f ') of the distribution, as a function of L in log-log 
scale for both type of disorder. Estimates for 9c are taken 
from Ea. (|22p and the straight lines indicating the asymptotic 
behavior have the same slope: — l/!^s ~ —0.8. The effective 
critical exponents calculated from Ea. (|16p are shown in the 
upper inset (for the mean) and in the lower inset (for the 
median). 



same for both types of disorder and can be well fitted as 
— 1/uu] ~ —0.8. As for the shift exponent in Sec lIV C"2l 
we have measured effective, size-dependent critical expo- 
nents, which are defined as: 



1 



Vw{L) 



= sinh 



Mc{2L) - Mc{L/2) 
2KejL) 



1 



(18) 



and plotted in the inset of FiglS] as a function of 1/L. 
Extrapolating the effective exponents yields 1/vw ^ 

0.805(10) and l/v^J'^ = 0.811(10) for the box-/i and the 
fixed-/i randomness, respectively. These values indeed 
agree within the error of the calculation, thus we can 
conclude that l/vyj ~ 0.808(10) and thus 



v^ = 1.24(2) 



(19) 



Comparing our estimates for the shift-exponent in 
Ea. p?|) with that of the width-exponent in Ea. ((T9| we 
notice that they agree within the error of the method, 
which corresponds to the rcnormalization group result 
for a classical conventional random critical pointi^. In 
order to make a direct check of the equivalence of the 
two exponents we have formed the ratio: 



a{L) = 



'c - 0c{L) 



(20) 



which should approach an i-indcpendcnt constant value 
at the 'true' critical point, 9c, provided Vs — ^w In 
FiglT] we have plotted the a{L) ratios as a function of 



InL using different input values for the critical point, 9c- 
As one can see in this figure the i-dependence of a{L) is 
very sensitive to the input value of 6c, both for box-/i and 
fixed-ft. randomness, but at its right value the a{L) ratios 
are approximately L independent. In this way we have 
demonstrated, that the infinite disorder fixed point of the 
2d RTFI model is characterized by one correlation-length 
exponent, which is given by: 



1.24(2) 



(21) 



Furthermore the ratio in Eq. (PH]) at the critical point has 
the universal value: a — 1.15(2), which does not depend 
on the form of the randomness. Finally, by this method 
we have obtained accurate estimates for the 'true' critical 
points, which are given by: 



6'('') = 1.6784(1) 
e^f'^ = -0.17034(2) 



(22) 



We have checked, that the values in Eg. (|22|) are consis- 
tent with other estimates, which can be obtained by ex- 
trapolating the 9c{L) data through Eq. p^ . but the error 
bars are smaller. For the box-/i disorder the known esti- 



3(&) 



(b) 



mates are 6*^ = 1.676(5), in Ref.[l8l and O)"' = 1.680(5), 
in Ref ■ |17| . which are consistent with that in Eg. d^^ . 
however the present value has much smaller uncertainty. 
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FIG. 6: (Color online) Standard deviation of the distribu- 
tion of the pseudo-critical points as a function of the size 
in log-log plot for box-/i (upper points) and fixed-/i (lower 
points) randomness. The dotted (blue) straight line has a 
slope: 1/vw = —0.808 corresponding to the estimated value. 
Inset: finite-size estimates for the exponent, l/u^, plotted as 
a function of l/L. 



V. SCALING AT THE CRITICAL POINT 

Having accurate estimates of the critical points for 
both types of randomness we are ready to study the crit- 
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FIG. 7: (Color online) The ratio in Ea. (piHl as a function of 
InL for different input values of the critical point, 6c- At the 
'true' critical point a(L) is approximately L independent. Up- 
per panel: box-/i randomness and 6c varies equidistantly be- 
tween 1.6786 and 1.6782 from up to down; lower panel: fixed- 
h randomness and dc varies equidistantly between —0.17030 
and —0.17038 from up to down. 



ical behavior of the system. In this respect we have con- 
centrated our effort at the critical point, where we have 
studied the distribution function of the magnetization, as 
well as that of the (log-)gaps and calculated critical expo- 
nents by finite-size scaling. At the critical point we have 
considered finite systems up to a linear length L = 2048 
and studied 4 x 10'' realizations for each sizes. 

These investigations are supplemented with numerical 
studies outside the critical point, both in the disordered 
and in the ordered phases. For the magnetization we have 
studied the scaling regime, defined as 8L^I^ = C(l). To 
obtain information about the dynamics of the system we 
have studied the scaling behavior of the excitation ener- 
gies both in the ordered and in the disordered Griffiths 
phases and in this investigation we have not restricted 
ourselves to the vicinity of the critical point. In the off- 
critical region we studied random samples with fixed-/i 
randomness up to a linear length L ~ 512 and for 10^ 
realizations. 



A. Magnetization 

1- Spin clusters 

During renormalization effective spin clusters are 
formed, which are than decimated at different energy 
scales. We illustrate the cluster structure of a given sam- 
ple at the critical point in the insets of FiglS] in which 
clusters having the same size are marked with the same 
grayscale (color) . Note that most of the clusters consist of 
only one site and the large clusters are generally geomct- 



rically disconnected. We have also analyzed the distribu- 
tion of the mass of the clusters, Pl{^), which follows the 
scaling form: L'^f P{fiL^'^f), df being the fractal dimen- 
sion defined in Eg. pTj) . According to scaling theorj*^ the 
distribution has a power-law tail, P(u) ^ u~^, with an 

exponent t = 1 H . The scaling function for L = 1024 

is shown in Fig. [8] in a log-log scale and indeed it has 
a linear dependence with a slope r = 2.9(1), which is 
consistent with our estimate for df in Eq. ([M)) . 




100 j^ 1000 



FIG. 8: (Color online) Distribution of the mass of the spin 
clusters in a log-log scale for the L = 1024 system at the 
critical point calculated from 40000 samples. The straight line 
indicating the asymptotic behavior has a slope — r = —2.9, see 
the text. In the two insets the cluster structure is illustrated 
for the L = 64 system with fixed-/i (upper inset) and box-/i 
(lower inset) randomness. The size of the clusters is increasing 
with the greyscale (color) as indicated at the bottom of the 
lower inset. The one-site clusters are white the points of the 
largest cluster are black. 

The magnetic properties of a given finite sample are 
related to the magnetic moment of some effective spin 
cluster, which appears at the last stages of the RG pro- 
cedure. In principle one can define different types of such 
spin clusters, i) The magnetization cluster has the small- 
est effective transverse field, thus decimated at the lowest 
energy scale. The moment of magnetization clusters are 
denoted by fl. By definition the smallest possible value 
of the magnetization cluster is flmin = 1, thus the corre- 
sponding minimal magnetization at the disordered phase 
is given by 1/L'*, i.e. it varies as a power-law of L. ii) 
The correlation cluster with a moment flcorr is present in 
the duplicated sample and involved in the replica corre- 
lation function. By definition a correlation cluster exists 
only for 9 < 9c{L), i.e. below the pseudo-critical point 
of the given sample. We have checked, that for large L if 
a correlation cluster exists, it is almost always the mag- 



netization cluster, iii) Finally we define also an energy 
cluster in a sample, which is the magnetization cluster 
for 9 > 9c{L), whereas for 9 < 9c{L) it is the cluster dec- 
imated before the correlation cluster. These clusters are 
involved in the low-temperature or the small longitudinal 
field properties of the system and will be considered in 
dynamical scaling in Sec lVBl 



2. Moment of magnetization clusters 

We have calculated the distribution functions of the 
moments of magnetization clusters for different lengths, 
Rl{JJ), which are shown in Fig[5] for both types of 
randomness. According to scaling theory Rl[j1) ~ 
L"^' ]k.{jiL~'^f) which is illustrated in the insets of FiglHl 
Up to a multiplicative constant the scaling functions, 
R{u}), are identical for the two different randomness 
and can be approximated with an exponential function: 
R{Ld) ^ exp(— w/o;*), uj* being some randomness depen- 
dent value. 



3. Fractal dimension and critical exponent 

In order to obtain an accurate estimate for the fractal 
dimension, df, and for the magnetization exponent, x, 
we have calculated average moments of the magnetiza- 
tion clusters, JI]^, which are plotted in the inset of Fig. 
[TUl as a function of L in a log-log scale. For both type 
of randomness the points are asymptotically on straight 
lines having the same slope, which is in agreement with 
the scaling relation in Eq. pH)) . We have calculated effec- 
tive, size-dependent fractal dimensions through: 



df{L) = sinh 



M2L Ml/2 



1 
ln2 



iTTT > (23) 



which are plotted in Fig. [TU] As seen in this figure the ef- 
fective fractal dimensions show no systematic trend with 
L and the df{L) values spread around the same mean 
value for both form of disorder. This mean value is taken 
to our estimate for the fractal dimension: 



df = 1.018(15) 



(24) 



and from Eg. ljlip we have for the magnetization expo- 
nent: 

a; = 0.982(15) . (25) 



4-. Scaling of the magnetization 

The magnetization is given by the asymptotic value of 
the correlation function in Eq.®, which definition can 
be extended in a finite system in terms of the replica 
correlation function as introduced in Sec lIVBl In a given 
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FIG. 9: (Color online) Distribution of the moment of mag- 
netization clusters for different lengths. Upper panel: box-/i 
randomness, lower panel; fixed-/i randomness. The scaled 
distributions are shown in the insets, where for the fractal 
dimension the estimate in Ea. (|24|) is used. 



0.95 



0.9 





f 






X X 


X 


+ 


X 






+ ^ X + V X 
_ _ _j^_ _ _^_ ^ _x- - - X- - - - - 4- - -X — -X 






+ 


+ + 








^ + 

+ 




b 




10= 
10^ 


f 

xx,---++ b 

X •' + 


- 










+ 






10^ 
10° 












10 1 00 1 000 








L 



10 



100. 



1000 



FIG. 10: (Color online) Finite-size estimates for the fractal 
dimension of the magnetization cluster in Eq. (|23|l as a func- 
tion of ln(L) for the Rxed-h (f) and the box-fe (b) randomness. 
The dotted horizontal line at df = 1.018 represents the mean 
value and our estimate. In the inset the average moment of 
the magnetization clusters are shown as a function of L in a 
log-log scale for both type of randomness. The dotted straight 
line has the slope df as extracted from the main figure. 
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sample of linear size L two spins at a distance r ^ L 
are correlated if both are in the same correlation cluster. 
Consequently (the average value of) the magnetization 



is given by: m{d, L) 



Ld 



-, where ^J.corr is the average 



value of the mass of the correlation cluster as defined in 
Sec IV A 1] The magnetization as a function of 5 is plotted 
in Figllllfor different finite systems. Using scaling theory 
the magnetization in the vicinity of the critical point is 
expected to behave as m{6,L) = L^^ 171(6 L^^"). To test 
this assumption in the inset of Fig ITT] we have plotted 
m{S, L)L^ as a function of 5L^/^ . Using the estimates for 
the critical exponents in Eas. ipT]) and ([25]) an excellent 
scaling collapse of the data is obtained. 



FIG. 11: (Color online) Magnetization as a function of the 
control-parameter in the vicinity of the critical point for dif- 
ferent finite systems. In the inset the scaled magnetization 
m{5, L)L^ is plotted as a function of SL^''^ . 



B. Dynamical scaling 

Here we study the properties of the low-energy exci- 
tations, which are responsible for the dynamical behav- 
ior of the system, such as the auto-correlation function 
or the low-temperature and small-field behavior of the 
susceptibility, specific heat, magnetization, etc. Such a 
low-energy collective excitation of a random sample of 
linear length L in the SDRG method is represented by 
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a large spin cluster which is formed at the last steps of 
the renormalization process and its excitation energy, e^, 
is measured by the effective transverse field of the clus- 
ter. These clusters are the so called energy clusters as 
defined in Sec IV A 11 According to the scaling relation 
in Eq.([5|) at the critical point it is convenient to use the 
log-variable: 7l = — In(eL). In the following we study 
the distribution of 7l at the critical point, as well as in 
disordered and ordered Griffiths phases, and investigate 
its scaling behavior with L. As far as dynamical prop- 
erties are considered the SDRG method gives asymptot- 
ically exact results also in the off-critical region, where 
the relaxation time is divergent^. This is not the case 
for static quantities due to the finite correlation length. 



1. Critical point 

At the critical point the energy clusters are decimated 
either at the last step of the renormalization process, 
which happens if the pseudo-critical point of the sam- 
ple, 9c{L), is smaller than 6c, or at the last but one step, 
if 9c{L) > 9c- The distribution of the log-excitation en- 
ergies at the critical point for different sizes are shown 
in Figd^l for the two type of randomness. As seen in 
this figure the distributions broaden with increasing L, 
which is a clear signal of infinite-disorder scaling. Re- 
ferring to Eq.([5]) we introduce the scaling combination: 
7 = (7i — jo)L~''f' , in terms of which the distributions 
collapse to the same curve provided the exponent is 
■0 sa 0.48. This is iUustratcd in the insets of FiglEl The 
constant term, 70, used in the fitting process is found to 
be 0(1) and has only a little influence on the value of the 
exponent ip. The scaling functions in the insets of FigfT2l 
have the same form for the two types of applied random- 
ness and have a heavier tail than in Id. The skewness 
values at i = 1024 are s — 0.82(1) for 2d to be com- 
pared with s — 0.64(1) in Id. 

In order to obtain more accurate estimate for the ex- 
ponent "0 we have considered the mean value, ^{L), as 
well as the standard deviation, Aj{L), of the distribu- 
tions. We have formed the differences, dl{L), which 
arc dl{L) = j{L) — j{L/2) for the mean and dl{L) = 
Aj{L) — A'-f{L/2) for the width of the distribution, re- 
spectively and which are plotted in the inset of Fig|T3las 
a function of L in a log-log plot. As seen in this figure 
the points for the two quantities and for the two-type 
of randomness are on parallel straight lines the slope of 
which is compatible with ip « 0.48. In the next step we 
have calculated finite-size effective exponents through the 
definition: 



^(L)^^ln 



dl(2L) 



dl{L) 



(26) 



which are plotted in FiglT3]as a function of L. As seen in 
this figure there is some systematic trend of the points up 
to L = 128— 192, but for larger L-s there are only fluctu- 
ations around a mean value what we use as an estimate 
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FIG. 12: (Color online) Distribution of the log- excitation en- 
ergies at the critical point for different finite systems. Upper 
panel: hox-h randomness; lower panel: fixed-/i randomness. 
In the insets scaling collapse of the distributions in terms of 
7 = (7 — 'yo)L~^ are shown, with ip = 0.48. The constant 
is 70 = ~l-5(4) for box-/i randomness and 70 = —0.4(1) for 
fixed-/i randomness, respectively. 



for the exponent: 



i) = 0.48(2) 



(27) 



This limiting value is the same for both type of ran- 
domness and describes well the scaling behavior both the 
mean value and the width of the distribution. 



2. Disordered Griffiths phase 

In the disordered Griffiths phase the energy clusters, 
which are related to the low-energy excitations in a large 
system are decimated at the last step of the renormal- 
ization process. The scaling behavior of the low-energy 
excitations here is different from that at the critical point. 
We illustrate the scaling behavior of the distribution of 
the log-excitation energies for different sizes in Fig|14lfor 
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FIG. 13: (Color online) Effective exponents ^{L) as a func- 
tion ln(I/) calculated from the mean or from the width of the 
distribution of the log-excitation energies for the two type 
of randomness. Inset: finite differences of the mean and the 
width of the log-excitation energy distribution as a function of 
L in a log- log scale, see text. The slope of the dotted straight 
line is given by V' = 0.48. 



that for large 7 the tail of \n.p(^) is linear and its slope, 
—d/zsi{L), can be used to obtam an independent esti- 
mate for the dynamical exponent. The measured values 
of d/ Zsi{L) are given in the caption of Fig lU and these 
are compatible vifith those calculated from the shift of the 
distributions. There are, however finite-size corrections 
for L < ^, where the correlation length at the studied 
value of 5 is of the order of ^ = 0(10^). 
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the fixed-/i randomness at a distance S = 9.64 x 10^'^ 
from the critical point. As seen in this figure the shape 
of the (logarithm of the) distribution functions is very 
similar for different L-s and the curves are merely shifted 
with ln(L). This follows from the assumption, that the 
typical value of the excitation energy in a system of size 
L scales as a power-law: 



eL-^L- 



(28) 



z being the dynamical exponent, which is a continuous 
function of the control parameter, S > 0^. Consequently 
the appropriate scaling combination is: 



1 = 1L~ z ln(L) - 70 



(29) 



in terms of which the distribution functions have a scaling 
collapse, provided the appropriate value of the dynami- 
cal exponent is used. This is illustrated in the inset of 
Fig |14l An estimate for the dynamical exponent from the 
shift of the distributions, Zsh{L), can be obtained from 
the optimal collapse of the data points for sizes L/2 and 
L. In FiglHwe have d/z^,,(512) = 2.3(2). If the low- 
energy excitations in the system are localized then the 
scaled distribution function is suggested^^ to be given 
by the Frechet distribution known from extreme value 
statistics^ in the form: 



7 - exp -7- + \n{d/z) . (30) 



lnp(7 - 70) 



Indeed the scaled distribution function in the inset of 
FigHHis well described by the function in Eg. ipn]) . where 
only one fitting parameter, 70 in Ea. (|29p is used. We note 



FIG. 14: (Color online) Distribution of the log- excitation en- 
ergies in the disordered Griffiths-phase (at 5 = 9.64 x 10^'^) 
in a log-lin scale for different sizes. The slopes of the 
straight lines indicating the tail of the curves are d/z = 
2.5(2), 2.1(2), 2.1(2) and 2.1(2), for L = 64,128,256 and 
512, respectively. In the inset the scaled distributions are 
shown in terms of the variable in Ea. (|29p with d/z = 2.3(2), 
which is well described by the Frechet distribution (full line) 
in Eg. pO)) . 



We have repeated the previous calculation for several 
values of 5 in the disordered Griffiths phase and calcu- 
lated estimates for the dynamical exponent both from the 
shift of the distributions and from the slope of the tail. 
These estimates obtained at different L-s are shown in a 
log-log plot in Fig. [15] One can notice that the finite-size 
corrections are stronger for small (5-s, where the correla- 
tion length is comparatively larger. According to scaling 
theory^^ the dynamical exponent for small 5 behaves as 



- « c5''^' 



(31) 



and divergent at (5 = 0. We have checked the rela- 
tion in Eq. ipT]) and indeed in Fig 115] one can identify 
an approximately linear part for 5 < 0.02 having a slope 
sa 0.6. This value is compatible with our previous esti- 
mates 1^-0 = 0.60(6) using results in Egs. ipT]) and (P7)) . 

Closing this section we note that the dynamical ex- 
ponent enters into the singularities of different physical 
quantities. For example at low-temperature the suscep- 
tibility behaves as: xi^) ^ y'/^-i and the specific heat 
has the form Cy ^ T'^/^ . More details about the scaling 
relations in the Griffiths phase can be found i n^^'^° . 
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FIG. 15: (Color online) Estimates for 2/z at different points 
of the disordered Griffiths phase in a log- log plot. The esti- 
mates are calculated either from the shift of the distributions 
or from their tail at various finite sizes. The straight line with 
a slope 0.6 indicates the asymptotic behavior as given in Eq. 
(|31|l . In the inset data for the largest system, L — 512, are 
compared with the scaling curve in Eq. (|31[l with c — 36.5 in 
linear plot. 



3. Ordered Griffiths-phase 

In the ordered phase of a large system with (5 < 
there is an infinite magnetization cluster, vifhich is deci- 
mated at the last step of the renormalization procedure. 
The low-energy excitations of the system here are re- 
lated to the energy cluster, which has a finite extent and 
which is decimated just before the magnetization clus- 
ter. The distribution of the log-excitation energies for 
different sizes are shown in Fig |16l for the fixed-/i ran- 
domness at 5 — —1.916 X 10"^. Comparing the distri- 
butions with that of the disordered Griffiths phase in 
FigHH one can notice that in both cases the distribu- 
tions are not broaden but shifted with an L-dependent 
amount. There are, however, several differences in the 
two figures. In the ordered Griffiths phase the finite-size 
effects are stronger, therefore we went up to L = 1024. 
More importantly, the shift of the distributions in the 
ordered Griffiths phase is slower than linear with ln(i). 
This is connected to the scaling result, that the typical 
value of the excitation energy, e^ , is related to the size of 
the system as: Ine^ ^ — In^' (i), thus the appropriate 
scaling combination is^ 



7 = 7i — A In ' (L 



70 



(32) 



which is to be compared with Eg. ip^ . Indeed using the 
variable in Eg. ([32]) the distributions show a scaling col- 
lapse, as illustrated in the inset of FiglTBl 

Also the shape of the scaled distributions are different 
in the two Griffiths phases. In the disordered Griffiths 
phase the distributions in the inset of Fig |14l approach 



FIG. 16: (Color online) Distribution of the log-excitation 
energies in the ordered Griffiths-phase (at 5 = — 1.916 x 10^^) 
in a log-lin scale for different sizes. In the inset the scaled 
distributions are shown in terms of the variable in Ea. (|32p . 



a linear asymptotics from above, on the contrary in the 
ordered Griffiths phase in the inset of Figlin] the points 
bend below a straight line. This is compatible with the 
scaling result, that asymptotically^: 



lnp(7) 



-l' 



(33) 



Our data in the inset of Figliniare still not in the asymp- 
totic regime but the tail of the distribution clearly de- 
creases faster than linear for the large sizes. 



VI. DISCUSSION 

The concept of infinite disorder fixed point has been 
introduced quite recently^ and its basic properties have 
been demonstrated in partially exact calculations in 
different one-dimensional random quantum^"— i^i^ and 
stochastic system9i2r22. In higher dimensions, in partic- 
ular in two dimensions the calculations are numerical and 
have only limited accuracyi^"— . In the present paper we 
have considered the prototypical model in 2d having an 
IDFP the random transverse-field Ising model and stud- 
ied its critical properties by a numerical implementation 
of the SDRG approach. As follows from the concept of 
IDFP our numerical results are expected to be asymp- 
totically exact. 

In our approach we have used a very efficient compu- 
tational algorithm which made us possible to treat sam- 
ples which are ten-times larger in linear size, compared 
with previous calculations. In this way we could reduce 
the effect of finite-size corrections and could also study 
off-critical properties, such as scaling functions and dy- 
namical scaling in the disordered and ordered Griffiths 
phases. 
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The main results of our investigations are the following. 
We have extended the finite-size scaling study for pseudo- 
critical points and from their distribution we have ob- 
tained precise estimate for the correlation length critical 
exponent, which has been shown to govern both the shift 
and the width of the distribution. Using different types of 
randomness we have demonstrated that the IDFP is uni- 
versal, the critical exponents as well as the critical scaling 
functions are independent of the randomness used in the 
calculation. We have also studied cross-over phenomena 
with respect of the linear size of the system as well as the 
type of randomness used in the calculation. The some- 
times large errors and deviations between the results of 
previous numerical studies are presumably due to cross- 
over effects. This can be seen, e.g. in FigfTSl where the 
estimates of the exponent tp have strong finite-size as well 
as randomness- type dependence. 

As a result of the larger samples and the good statis- 
tics of the numerical data we have obtained accurate esti- 
mates for the critical exponents and studied - at the first 
time - the behavior of scaling functions, both at the criti- 
cal point and in the finite-size scaling limit in the vicinity 
of the critical point. We have also extended our investi- 
gations to the disordered and the ordered Griffiths phases 
and have checked various predictions of phenomenologi- 
cal scaling theory. 

Comparing the critical behavior in Id to that in 2d we 
have qualitatively similar results, but there are also im- 
portant differences. First of all the actual values of the 
critical exponents as well as the form of the scaling func- 
tions are different. Scaling in the ordered Griffiths-phase, 
which involves powers of d in Eqs. p2p and p3p however, 
is also qualitatively different. One particular feature of 
the model in Id is self-duality, which could be the reason 
why in Id the distribution of the pseudo-critical points 
involves different shift and mean exponents^. 

Our calculations can be extended to several directions. 
One possibility is to study the entanglement properties 
of the model^ which can be well performed within the 
frame of the SDRG approach^. At present conflict- 
ing theoretical predictions are available about the finite- 
size dependence of the critical entanglement cntropyi^iii, 
which could be possibly clarified by using larger systems 
in the calculation. A second promising direction of appli- 
cation of our approach is to consider higher dimensional 
systems. At present even in three dimensions only the 
existence of infinite disorder scaling is dcmonstratedi^, 
but no estimates are known about the critical exponents. 
In four and higher dimensions no numerical studies have 
been performed so far. Studies in these directions are in 
progress. 
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Appendix A: Filtering out irrelevant bonds 

In the SDRG procedure if the maximum rule is applied 
several bonds are irrelevant, which means that they can 
be deleted without modifying the final result of the renor- 
malization. 



Condition for a bond to be irrelevant 



Let us consider a bond with a log-coupling k^- — In Jij 
between sites i and j and consider such nearest neighbor 
points, one of those is denoted by k in Fig. 1171 which 
has bonds both to i and j of strength Kik and Kjk-, re- 
spectively. The triangle {i,j, k) is called the 'majorating 
triangle' of the (i, j) bond if Ky is the smallest bond in the 
triangle, and it is also smaller, than the potentially gen- 
erated new bond: k'^ = Kik + Hjk — ^fc, where 9^ = In/ifc 
is the log-transverse field at site k. If such a majorating 
triangle exists then the bond is irrelevant. 



2. Proof of the filtering criterion 

For the proof it is enough to consider such a decimation 
step during which the majorating triangle collapses and 
we follow also the evolution of the parameters in another 
triangle, {i,j,l), which is not necessarily a majorating 
one. The maximal log-parameter in the four-site system 
in Fig. [T7]is denoted by uj. We note that during renor- 
malization the log-transverse fields can not increase, thus 
the change of their values are unimportant in the proof. 



a. Bond decimation 

The (i,j) bond is not the largest one, since it has a 
majorating triangle, so that for the position of the largest 
bond we have three different cases, i) Largest bond in the 
majorating triangle: oj = Kik (or Kjk). The spins i and k 
fuse into a new effective spin, which is connected to the 
spin j with a coupling of k' = max (k^ , Kjk) = Kjk, thus, 
Kij simply disappears, ii) Largest bond in the another 
triangle: uj = kji (or ku). The j and I spins fuse into 
a new effective spin. If Kij < ku, than the value of Ky 
simply disappears, otherwise its value does not change. 
Similarly, the value of Kik does not change, while Kjk is 
replaced by k' j. = max {Kjk, Kki) > Kjk- Taken all round, 
the new (i,j,k) triangle is a majorating triangle of the 
(i,j) bond, iii) Largest bond between the two triangles: 
ll) — Kki ■ The spins k and I fuse and there is one triangle 
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FIG. 17: Illustration of the neighborhood of an irrelevant 
bond {i,j) having a majorating triangle {i,j,k). The values 
of the log-couplings and the log-transverse fields are also in- 
dicated, see the text. 



left. Here the value of k^ does not change and the new 
triangle is a majorating one of the (i, j) bond. 



Site decimation 



We should consider two different cases: \) uj ~ Ok'- 
K^ = Kife -|- Kjk ^ Ok > Kij due to our assumptions, thus 
Kij disappears without affecting the results, ii) lo = 9i (or 
9j analogously): k'-j. = k^- -|- Ki^ — 9i < Kij < Kjk, thus 
this coupling is unaffected by the k^- coupling. How- 
ever, in principle the other newly generated bonds are 
not irrelevant, e.g. «;'; < Kji is not always fulfilled. Due 



to this, the effect of the ma j orated (?, j) edge will not 
certainly disappear after one single decimation step, in 
which it is involved. Moreover, a lot of new couplings can 
be generated occasionally, if the bond lies in a 'dangerous 
position'. 



c. Completing the proof 



Here we show, that the newly generated coupling in 
the {l,j,k) triangle, k';, is always majorated by this tri- 
angle. The two neighboring couplings of k'; are Kjk (due 
to the fact, that k'j, < Kjk), and max (Kki , K'ki) , where 



^kl 



Hik + l^il 



is smaller than these neighbor- 



ing couplings, namely k'j < Kij < Kjk and kJ, 



^kl 



^Jl 



K-ik — Hij > 0. From the latter follows, that k'-; < k'^j, 
which can not be greater, than max {Kki, kJ,;) correspond- 
ing to the new value of the coupling between the k and I 
spins. Now we see, that the k'-^ coupling is majorated by 
the (l,j, k) triangle, if k'-^ is smaller, than the generated 
coupling k'ji = max {km, K'j,i) + Kjk — 9k, which is obtained 
by decimating the transverse field at k. Let us consider 
their difference: k'L — k'-^ > k'j,j + Kjk — 9k — Kij — ku + 9i = 
Kik + Kjk — 9k — Kij = Kij — Kij > 0. Thus the k'j coupling 
is always majorated by the {l,j, k) triangle. 

Let us summarize our findings. If the bond {i,j) has 
a majorating triangle, it can not be decimated directly 
Decimating in its neighborhood either this bond disap- 
pears or new couplings are generated, the value of which 
involves Kij. These new couplings, however, are always 
majorated by a triangle, thus they are never decimated 
during the renormalization process. Consequently the 
value of Kij does not influence the result of the renormal- 
ization. 
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